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Elastic properties of the filled skutterudite compound SmOs4Sbi2 have been investigated by ultra- 
sonic measurements. The elastic constant Cn(&>) shows two ultrasonic dispersions at ~15 K and ~53 
K for frequencies ll> between 33 and 316 MHz, which follow a Debye-type formula with Arrhenius-type 
temperature-dependent relaxation times, and remain unchanged even with applied magnetic fields up to 
10 T. The corresponding activation energies were estimated to be E2 = 105 K and E l = 409 K, respec- 
tively. The latter, E\ , is the highest value reported so far in the Sb-based filled skutterudites. The presence 
of magnetically robust ultrasonic dispersions in SmC^Sbn implies a possibility that an emergence of a 
magnetically insensitive heavy fermion state in this system is associated with a novel local charge degree 
of freedom which causes the ultrasonic dispersion. 
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Intermetallic compounds with cage like structures have 
been studied both experimentally and theoretically for several 
decades. 1 ' 2 ^ Recently, this area of research has received in- 
tense interest, since a variety of exotic physical properties are 
observed in these types of compounds; e.g., heavy-electron 
states of heavier rare-earth (Pr, Nd, and Sm) and La based 
systems, multipole ordering, unconventional superconductiv- 
ity, high thermoelectric figures of merit, and anharmonic vi- 
brations of the guest ions. ' It has been suggested that these 
exotic phenomena may originate from the characteristic cage 
structure in these systems. A well-known system of inter- 
metallic system, with a cage like structure, that exhibits a va- 
riety of correlated electron phenomena, are the filled skutteru- 
dite. 3_6) For instance, an unconventional heavy fermion (HF) 
state is observed for SmOs4Sbi2, T) where the electronic spe- 
cific heat coefficient -y=820 mJ/K 2 mol is insensitive to an ap- 
plied magnetic field H || [100] up to 8 T. 8) This result is in con- 
trast to the expected suppression of y with magnetic field for 
the conventional (magnetic) Kondo effect, where the heavy 
fermion ground state is derived from spin degrees of freedom, 
as is often observed for Ce- and Yb-based compounds. In con- 
trast, the exotic HF state in SmOs 4 Sbi2 might be explained 
by a novel many-body effect, such as the multichannel Kondo 
effect, which is derived from multipole degrees of freedom 
and/or local charge fluctuations originated from off-center de- 
grees of freedom of guest ions. 9 ' 10) 

By means of ultrasonic measurements, a frequency depen- 
dent upturn of the elastic constant with ultrasonic attenuation 
peak, called an ultrasonic dispersion (UD), has been observed 
in the cage compounds, such as R^P^Gtf, (R = La, Ce, Pr) 
and /?Os 4 Sbi2 (R = La, Pr, Nd). 11_13) Moreover, an unusual 
decreasing of the elastic constant has also been observed at 
very low temperatures in the La based compounds. Due to the 
co-presence of UD and low temperature elastic softening in 
non-4/ La based compounds, the origins of these phenomena 



are expected to share a common root cause, originated with a 
novel local charge degrees of freedom. As a candidate for the 
origin, a model based on an off-center ionic configurations in 
a multi-well potential has been proposed. Here, the UD can be 
understood as a thermally activated off-center motion of guest 
ion, which was defined as e(off-center) rattlingf in our previ- 
ous papers. 13 ' I4) On the contrary, recent neutron scattering ex- 
periment on PrOs 4 Sbi2 has revealed that the Pr-nuclei density 
distribution with an accuracy of 0.1 A is on-center-like with- 
out anisotropy at 8 K while it is off-center-like with strong an- 
harmonicity at room temperature. 15 ' Therefore, it is still con- 
troversial whether the off-center configuration of guest ion is 
appropriate for understanding the above elastic anomalies. 

Here, we note that the word 'rattling' has originally been 
used for the localized thermal vibration of the guest atom in 
an over-sized host cage with large amplitude. 6 * In the rattling 
systems, a low-lying optical phonon excitation of ftai ~ meV 
has been commonly observed in several physical quantities, 
e.g., Raman scattering, inelastic x-ray or neutron scattering 
and Einstein temperatures in the specific heat. 16-19 ' On the 
other hand, since the ultrasonic measurement with frequen- 
cies of ~ MHz can monitor relatively lower energy excita- 
tions of fid) ~ 10~ 5 meV, a causal linkage between the op- 
tical phonon and UD has not been confirmed thus far and 
still an open question. In order to shed light on the nature 
of the rattling phenomenon and UD, a systematic change of 
the characteristic parameters obtained from the UD and other 
physical properties among the Sb-based filled skutterudites 
will be reviewed in the present paper. In addition, although 
ultrasonic measurements of SmOs4Sbi2 have already been re- 
ported, 20 ' there is no mention of the ultrasonic frequency de- 
pendence. The present study is focused on investigating the 
UD in SmOs4Sb 12 in order to examine the magnetically robust 
novel local charge degree of freedom, which might be a possi- 
ble candidate for the origin of exotic HF state in SmOs4Sbi2. 
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Fig. 1. Elastic consants Cu (176 MHz) and C 44 (165 MHz) of SmOs 4 Sb 12 
as a function of temperature T, displayed as a relative change for compar- 
ison. Inset: the green and red curves, displayed on a logJ scale, represent 
the fits of the quadrupole susceptibility by using the Tj and Fg, ground 
states, respectively. The backgrounds of the elastic constant C^ 4 curves for 
the fits are displayed as the same colored-dotted lines in the main panel. 

Single crystals of SmOs4Sbi2 were grown using a molten 
flux growth method. The high quality of the single crystals 
used in the present study were confirmed by energy disper- 
sive x-ray spectroscopy, powder x-ray diffraction (for single 
crystals picked out of same sample batch), the large residual 
resistivity ratio (RRR (300 K/ 2 K)) of 16.5, and an observa- 
tion of acoustic de Haas-van Alphen signals in the elastic con- 
stant C\ i (not shown). A cubic shaped specimen, with a length 
of 0.775 mm, was used for the ultrasonic measurements. A 
change in the sound velocity was detected by the conven- 
tional phase comparator method with pulsed ultrasound gen- 
erated by LiNb03 transducers. Ultrasonic attenuation mea- 
surements were performed by recording the in-phase-(sine) 
and quadruture-(cosine) signals at fixed frequencies and the 
phase shift. 21 ' An Oxford Instruments Heliox-TL 3 He refrig- 
erator with a superconducting magnet was used for measure- 
ments at T > 300 mK and H < 10 T. 

A comparison of the elastic constants C\\ and C44, as a 
function of temperature, is displayed as the relative change in 
ACij/Cjj in Fig. 1. Note that C u = C B + f C " 2 C ' 2 consists of 
ri-symmetry bulk modulus Cb = c " + 3 2c ' 2 and r23-symmetry 
c "~ c ' 2 mode while C44 mode associates with ^-symmetry. 
Cu shows two upturns that are frequency-dependent and 
magnetic field independent (details will be described later in 
Figs 2 and 3). The temperatures for the appearance of these 
two UDs are marked by vertical dotted lines in Fig. 1. In con- 
trast, C44 shows neither an upturn nor a frequency dependence 
at these temperatures, except for the elastic softening due to 
the crystalline electric field (CEF) effect from ~ 100 K to 10 
K and a step-like decrease at Tq = 2.5 K, below which a weak 
ferromagnetic moment develops for SmC^Sb^. 7 ' 

Similar contrasting behavior between C\\ and C44 was also 
observed in the /?Os4Sb;2 (R = La-Nd) systems. These ultra- 
sonic mode-dependence suggest that the UD would be derived 
from a symmetrical off-center charge fluctuation of the guest 
ion, e.g., r23-symmetry off-center mode for/?Os4Sbi2 and T$- 
symmetry for/?3Pd2oGe6. I2) We note that our results, the soft- 
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Fig. 2. (a) Relative change in the elastic constant ACn/Cn as a function 
of temperature T with frequencies varying from 33 to 316 MHz. (b) The- 
oretical fits to the frequency-dependent elastic constant. Arrows indicate 
the temperatures satisfying the resonant condition cot ~ 1, (c) Ultrasonic 
attenuation coefficient au as a function of temperature. Solid lines repre- 
sent phenomenological calculations using the same parameters as used in 
the calculations in (b). 

ening of C44 in particular, are completely different from pre- 
vious work by Nakanishi et al. 2ff> The major difference could 
simply come from the fact that higher frequency and shorter 
pulses, i.e., more focused ultrasonic waves, were used in the 
present measurements, which help to minimize the irregular 
reflection of ultrasonic waves from the crack or cavity on sur- 
face or inside the crystal. A good example of this effect can 
be seen in a deviation of the lowest frequency (33 MHz) data 
of C\\ in Fig. 2 (a). 

The softening of C44 can be analyzed in terms of 
quadrupole susceptibility based on the CEF effect. 22 ' Thus far, 
two different CEF energy level schemes have been proposed 
for the Sm 3+ ion (J = 5/2) in SmOs4Sbi2 in Oh symmetry: 
i.e., F 7 (0 )-r 8 (38 ) and T 8 (0 ) - F 7 (19 ). 7 ' 23) The inset of Fig. 
1 shows theoretical fits to the experimental data below 200 K 
for both CEF schemes, both of which reproduce the softening 
of C44 down to ~ 10 K. The leveling off of C44 below ~ 10 K 
is described by Van Vleck-type behavior, derived from a r 7 - 
doublet ground state using fitting parameters including a cou- 
pling constant for a quadrupole-strain interaction \g\ = 258 K 
and for an inter-site quadrupole interaction g' = -0.79 K. In 
contrast, the scheme with a F^-quartet ground state gives a fit 
which deviates from the C44 data due to a Curie-type contri- 
bution from the Fg-quartet ground state, that persists down to 
Tc- We should, however, note that the large value of the low 
temperature specific heat indicates that complex many-body 
physics dominates the low temperature state in SmOs4Sbi2- 8) 
Therefore, it is likely that the deviation of the C44 data from 
the fit, i.e., the shoulder for T < 10 K, should probably be ex- 
plained by a modified quadrupole susceptibility picture that 
includes a crossover from a Kondo singlet state (e.g., as de- 
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scribed for CePc^A^) ' ' or some other exotic Kondo state 
(e.g., the multichannel Kondo effect). 10,26) Thus, the possibil- 
ity of the iVquartet ground state still remains. If we skip the 
adjustment of the many -body effect, the fitting parameters for 
the r§ model assume relatively larger values: \g\ = 369 K and 
g' = — 1 .76 K. In the present paper, we will focus on the ul- 
trasonic dispersions and leave the discussion of the correlated 
electron physics on the CEF effects for a later publication. 

Figures 2 (a) and (c) show the relative change in the elas- 
tic constant ACn/Cn and the ultrasonic attenuation coeffi- 
cient a ii as a function of temperature with fixed frequencies 
varying from 33 to 316 MHz. The two upturns in C\\, indi- 
cated by the arrows in Fig. 2(a), are found to shift to higher 
temperatures with increasing frequency. Maxima in an are 
also observed at ~ 53 K and ~ 15 K, which show the same 
frequency dependence as the upturns in C\\. The frequency 
dependent elastic constant C(co), to is ultrasonic frequency, 
and the attenuation coefficient a(io) are well described phe- 
nomenologically, as shown in Fig. 2(b), by a Debye-type dis- 
persion given by, 14) C(w) = C (oo ) - Xl;=i,2 and a(a>) = 

ao - 2, =1,2 2~p~ i+J^ T 2 ■ where an Arrhenius-type temperature 
dependence of the relaxation time, r, = T^(i)exp(Ej/k B T), 
is assumed. Here, i = 1 and 2 are the indices for analyzing 
two separated sets of UDs at ~ 53 K and ~ 15 K, respec- 
tively. Here, AC, in the above equations is the total change 
for each upturn between the high-frequency limit Coo and the 
low-frequency limit Co of the elastic constant, which are esti- 
mated from the extrapolation of the temperature dependence 
of Cn over 20 K, without the UDs and possible CEF effect (as 
seen below 10 K in Fig, 2(a)), as displayed as dotted curves in 
Fig. 2(b). Two parameter sets of activation energy E, and at- 
tempt time ro,(i) were obtained from the phenomenologicalfits 
for the relaxation cot ~ 1 of Fig. 2(b) and 2(c), even though 
the ultrasonic attenuation peaks are often experimentally ob- 
served sharper than the phenomenological fits. 14 ' The param- 
eters are summarized in Table 1 with the related compounds' 
parameters of UD and also other physical properties, includ- 
ing the lattice constant a, the Debye and Einstein temperatures 
<9d and 6>e, and the Sommerfeld coefficient y. 

Figure 3 shows Cn vs T at a frequency of 105 MHz and 
for various magnetic fields between and 10 T. The UDs are 
unaffected by magnetic fields up to 10 T except for the field 
dependent CEF effect that appears in Cn below 15 K, as can 
be seen in the enlarged scale in the inset of Fig. 3. Such insen- 
sitivity in the UD to magnetic fields is commonly observed in 
the /?Os4Sbi2 (R — La-Nd) system, providing strong evidence 
that the UD is not magnetic but rather electric, i.e., the novel 
local charge degrees of freedom, in origin. 

Figures 4 (a) and (b) show a systematic change of the UDs 
and an Arrhenius plot of the relaxation time r, vs l/T for 
/?Os4Sbi2 (R = La-Sm), respectively. The double UDs so far 
only have been observed in NdOs4Sbi2 (4/ 3 ) and SmOs4Sbi2 
(4/ 4 ), while LaOs 4 Sbi2 (4/°) and PrOs 4 Sb 12 (4/ 2 ) show sin- 
gle UD. It has been pointed out that the observations of dou- 
ble UDs are comparable to the 123 Sb-NQR measurements on 
/?Os4Sbi2 (R = La-Sm), which show double or triple peaks 
in the temperature dependence of I/T2 for R = Pr-Sm due 
to some local charge fluctuation with a correlation time of 
~ 10~ 6 sec, while LaOs4Sbi2 shows only a single peak. 27 ' 







0T 

\ 21 
\ 57 


SmOs 4 Sb 12 - 






\ 10 T 


k || u || [100] || H " 
V 105 MHz 


(80 K) 


■ r ^z 


\ f 

\^ 0.1% 


" O 


















1 T(K) 10 













~20 4cT~ 60 80 
T(K) 



Fig. 3. Temperature dependence of the relative change of C\\ for various 
magnetic fields H along the [001] direction. The inset displays an enlarged 
view below 20 K on a log7" scale. 
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Fig. 4. (a) Comparison of the UDs that appear in elastic constant Cn of 
/?Os4Sbi2 (R = La-Sm) at several frequencies. Lower arrowheads with 
numbers 1 and 2 indicate the relaxation point air, ~ 1 for 1 = 1 and 2, 
respectively. The displayed data have been shifted to eliminate overlap- 
ping with each other and the SmOs4Sbi2 data are magnified three times 
for the ACn/Cii-axis. (b) Arrhenius plots of the characteristic parameters 
of the UDs (t (,> E t ) for i?Os 4 Sbi 2 (R = La-Sm). 

When the rare-earth ion changes from R — La to Sm, not 
only the guest ion's mass but also the amount of free space for 
the Einstein oscillator in the atomic cage increases, due to the 
lanthanide contraction, i.e., the ionic radius of the guest ion 
decreases while the cage radius remains nearly constant for 
7?Os4Sbi2. 6) We speculate that such changes in the guest ion's 
conditions give rise to additional anharmonicity of the poten- 
tial, e.g., flat base or double well, for the local ionic configu- 
ration, and result the emergence of additional UD and NQR 
peaks in R — Nd and Sm. 

In Fig. 4 (b), the slope of the line is equivalent to the acti- 
vation energy and the markers indicate frequencies and tem- 
peratures where the relaxation of UD meets a resonant con- 
dition u>Tj ~ 1. We find that the activation energy Ei and 
Sommerfeld coefficient y increase and ro,(i) decreases from R 
= La to Sm in the 7?Os4Sbi2 series as summarized in Table 1. 
The increase y means that a characteristic temperature of the 
electronic correlation, such as Kondo temperature, will ex- 
hibit opposite (decrease) tendency in general. Therefore, the 
similar increase tendency of E\, E2 and y, found in the present 
study, implies that the activation energies are not directly con- 
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Table I. Comparison of characteristic parameters of UDs and rattling, and some physical properties of the related filled skutterudites. 6,12 ' 19 ' 28 ' 





Lattice Constant 


Debye Temp. 


Einstein Temp. 


Sommerfeld Coef. 


Activation Energy 


Attempt Time 








a (A) 


@D (K) 


@E (K) 


y (mj mol^'Kr 2 ) 


Ei (K) 


T0,(D (P s ) 


E 2 (K) 


T0,(2) (ps) 


LaOs4Sbj2 


9.3081 


270 


60.5 


36-56 


127 


50 






CeOs 4 Sbi2 


9.3011 


304 




92-180 










PrOs 4 Sb 12 


9.3031 


165-320 




310-750 


225 


31 






NdOs 4 Sbi2 


9.2989 


255 


39 


520 


337 


7.5 


67 


51 


SmOs4Sbi2 


9.3009 


294 


40.1 


820 


409 


4.4 


105 


10 


LaRu 4 Sbi2 


9.2700 


275 


72.8 


47 










LaFe 4 Sbi2 


9.1395 


314 


87.6 


122 


300 


31 






PrFe 4 Sb 12 


9.1351 






-1000 


360 


3 







nected to the characteristic temperature of the electron mass- 
enhancement mechanism in SmOs4Sbi2- 

Next we compare the Einstein temperature with activation 
energies. If we pick out the lower activation energy of each 
compound for comparison, i.e. E\ = 127 K for La, Ei = 67 K 
for Nd and E 2 = 105 K for Sm, they seem to be qualitatively 
comparable to the systematic change of Einstein temperature 
6> £ ; 60.5 K for La, 39 K for Nd and 40.1 K for Sm in Table 
1 . This fact implies that the optical phonon and UD could be 
linked with unknown pre-factor of ~ 2. This linkage might be 
due to a novel energy dissipation mechanism via conduction 
electron-phonon interaction, as previously proposed by Hat- 
tori and Miyake. 10) If there exists such conduction electron- 
phonon interaction, it can not be ruled out that the origin of 
activation energies is indirectly involved in the effective mass- 
enhancement mechanism, which might be mediated by the 
interaction between conduction electron and phonon in the 
present system. 

Finally, we discuss a possible co-presence of Curie-type 
softening and UD in SmOs4Sbi2- The low temperature soft- 
ening as observed in LaOs4Sbi2 could, however, not be found 
in the present study of SmOs4Sbi2 because the drastic change 
of elastic constants due to ferromagnetic ordering at 2.5 K 
and the CEF effect under magnetic fields hides such features. 
In order to elucidate the possible ground states of the local 
charge degrees of freedom and extract the CEF level scheme 
in the present compound, further measurements, such as elas- 
tic constant (C\\-C\i)l2 and C44 measurements in a magnetic 
field, are needed. In the present discussions, the Sm-ion's va- 
lence fluctuations 2 ^ are not taken into account. We can, how- 
ever, conclude that the magnetically robust HF compound 
SmOs4Sbi2 at least has the magnetic field-independent relax- 
ation of local charge degree of freedom which cause ultra- 
sonic dispersions and will be a key to understanding the novel 
many-body effect in the present system. 

This work was supported by Grant-in-Aid for Scien- 
tific Research on Innovative Areas 'Heavy Electrons' (No. 
21102501) and for Specially Promoted Research (No. 
18002008) 'Strongly correlated quantum phases associated 
with charge fluctuations' of the MEXT, Japan. Research at 
UCSD was supported by the U.S. Department of Energy un- 
der Grant No. DE-FG02-04ER46105 and the National Sci- 
ence Foundation under Grant No. DMR 0802478 for single 
crystal growth and characterization. One of the authors (T.Y.) 
was supported by Hokkaido Univ. Leader Development Sys- 
tem in the Basic Interdisciplinary Research Areas and (R- 



C.H.) was supported by Research Corporation under CCSA 
Grant No. 7669. 



1) V. Keppens, D. Mandrus, B. C. Sales, B. C. Chakoumakos, P. Dai, R. 
Coldea, M. B. Maple, D. A. Gajewski, E.J. Freeman, and S. Bennington: 
Nature (London) 395 (1998) 876. 

2) C. C. Yu and P. W. Anderson: Phys. Rev. B 29 (1984) 6165. 

3) B. C. Sales: in Handbook on the Physics and Chemistry of Rare Earths, 
ed. K. A. Gshchneider, Jr., J.-C. Bunzli, and V. K. Pecharsky (north- 
Holland, Amsterdam, 2003) Vol. 33, p. 1. 

4) B. C. Chakoumakos, B. C. Sales, D. Mandrus and V. Keppens: Acta Crys- 
tal. B 55 (1999) 341. 

5) H. Sato: J. Magn. Magn. Mater. 310 (2007) 188. 

6) D. J. Braun and W. Jeitschko: J. Less Common Met. 72 (1980) 147. 

7) W. M. Yuhasz, N. A. Frederick, P.-C. Ho, N. P. Butch, B. J. Taylor, T. 
A. Sayles, M. B. Maple, J. B. Betts and A. H. Lacerda, P. Rogl and G. 
Giester: Phys. Rev. B 71 (2005) 104402. 

8) S. Sanada, Y. Aoki, H. Aoki, A. Tsuchiya, D. Kikuchi, H. Sugawara and 

H. Sato: J. Phys. Soc. Jpn. 74 (2005) 246. 

9) T. Hotta: J. Phys. Conf. Ser. 150 (2009) 042061. 

10) K. Hattori, Y. Hirayama and K. Miyake: J. Phys. Soc. Jpn. 74 (2005) 
3306. 

11) T. Goto, Y. Nemoto, T. Yamaguchi, M. Akatsu, T. Yanagisawa, O. Suzuki 
and H. Kitazawa: Phys. Rev. B 70, 184126 (2004) . 

12) Y. Nemoto, T. Yanagisawa, Y. Yasumoto, H. Kobayashi, A. Yamaguchi, 
S. Tsuduku, T. Goto, N. Takeda, A. Ochiai, H. Sugawara, H. Sato and H. 
Kitazawa: J. Phys. Soc. Jpn. 77, Suppl. A (2008) 153. 

13) Y. Nemoto, T. Yamaguchi, T. Horino, M. Akatsu, T. Yanagisawa, T. Goto, 
O. Suzuki, A. Donni and T. Komatsubara: Phys. Rev. B 68 (2003) 1 84109. 

14) T. Yanagisawa, Pei-Chun Ho, William M. Yuhasz, M. Brian Maple, Y. 
Yasumoto, H. Watanabe, Y. Nemoto and T. Goto: J. Phys. Soc. Jpn. 77 
(2008) 074607. 

15) K. Kaneko, N. Metoki, H. Kimura, Y. Noda, T. D. Matsuda and M. Kohgi: 
J. Phys. Soc. Jpn. 78 (2009) 074710. 

16) N. Ogita, T. Kondo, T. Hasegawa, Y. Takasu, M. Udagawa, N. Takeda, K. 
Ishikawa, H. Sugawara, D. Kikuchi, H. Sato, C. Sekine and I. Shirotani: 
PhysicaB 383 (2006) 128. 

17) S. Tsutsui, H. Kobayashi, J. P. Sutter, H. Uchiyama, A. Q. R. Baron, Y. 
Yoda, D. Kikuchi, H. Sugawara, C. Sekine, I. Shirotani, A. Ochiai and H. 
Sato: J. Phys. Soc. Jpn. 77 Suppl. A (2008) 257. 

18) K. Iwasa, Y. Mori, L. Hao, Y. Murakami, M. Kohgi, H. Sugawara and H. 
Sato: J. Phys. Conf. Ser. 92 012122. 

19) K. Matsuhira, C. Sekine, M. Wakeshima, Y. Hinatsu, T. Namiki, K. 
Takeda, I. Shirotani, H. Sugawara, D. Kikuchi and H. Sato: J. Phys. Soc. 
Jpn. 78 (2009) 124601. 

20) Y. Nakanishi, T. Tanizawa, T. Fujino, P. Sun, M. Nakamura, H. Sugawara, 
D. Kikuchi, H. Sato and M. Yoshizawa: J. Phys. : Conf. Ser. 51 (2006) 25 1 . 

21) B. Wolf, B. Luthi, S. Schmidt, H. Schwenk, M. Sieling, S. Zherlitsyn and 

I. Kouroudis: Physica B 294-295 (2001) 612. 

22) T. Goto and B. Luthi: Adv. Phys. 52 (2003) 67. 

23) Y. Aoki, S. Sanada, H. Aoki, D. Kikuchi, H. Sugawara and H. Sato: Phys- 
ica B 378-380 (2006) 54. 

24) Y. Shimizu and O. Sakai: J. Phys. Soc. Jpn. 65 (1996) 2632. 

25) Y. Nemoto, T. Goto, S. Nakamura and H. Kitazawa: J. Phys. Soc. Jpn. 65 



J. Phys. Soc. Jpn. 



Letter 



T. Yanagisawa et at 



5 



(1996) 2571. 

26) D. L. Cox and A. Zawadowski: Adv. Phys. 47 (1998) 599. 

27) H. Kotegawa, K. Tabira, Y. Irie, H. Hidaka, T. C. Kobayashi, D. Kikuchi, 
K. Tanaka, S. Tatsuoka, H. Sugawara and H. Sato: J. Phys. Soc. Jpn. 77 
Suppl. A (2007) 90. 



28) I. Ishii, H. Higaki, S. Morita, I. Mori, H. Sugawara, M. Yoshizawa and 
T. Suzuki: Physica B 403 (2008) 887. 

29) M. Mizumaki, S. Tsutsui, H. Tanida, T. Uruga, D. Kikuchi, H. Sugawara 
and H. Sato: J. Phys. Soc. Jpn. 76 (2007) 053706. 



